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The  work  described  herein  was  sponsored  by  Headquarters,  U.S.  Army  Corps 
of  Engineers  (HQUSACE),  as  part  of  the  Civil  Woiks  Research  and  Development 
Program  on  Structural  Engineering  (CWR&D).  The  research  was  conducted 
under  Civil  Works  Work  Unit  31589,  the  “Computer-Aided  Structural  Engineer¬ 
ing  (CASE)”  Project  for  which  Dr.  Stanley  C.  Woodson,  Geosciences  and 
Structures  Laboratory,  Vicksburg,  MS,  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC),  is  Problem  Area  Leader  and  Mr.  H.  Wayne  Jones, 
Computer-Aided  Engineering  Division  (CAED),  Information  Technology 
Laboratory  (TIL),  Vicksburg,  MS,  ERDC,  is  the  Principal  Investigator.  The 
HQUSACE  Technical  Monitor  is  Mr.  Jerry  Foster,  CECW-ED. 

The  SOILSTRUCT-ALPHA  for  Personal  Computers  analysis  package  com¬ 
prises  three  separate  stages:  the  Visual  Modeler,  the  SOILSTRUCT-ALPHA 
Finite  Element  Analysis  Program,  and  the  Visual  Postprocessor.  This  report. 
Report  1  of  a  series,  describes  how  the  SOILSTRUCT-ALPHA  Visual  Modeler 
(i.e.,  a  preprocessor)  is  used,  applying  it  to  the  suggested  work  flow  of  a  typical 
project  as  an  example.  This  report  is  intended  to  serve  as  a  user’s  guide  and  not 
as  a  reference  manual  defining  each  specific  command.  Report  2  will  describe  the 
Finite  Element  Analysis  Program,  and  Report  3,  the  Visual  Postprocessor. 

The  work  was  performed  at  ITL  by  Mr.  Barry  C.  White,  Computer  Systems 
Technology,  Vicksburg,  MS,  and  Dr.  Robert  M.  Ebeling,  CAED.  Dr.  Ebeling 
was  the  author  of  the  scope  of  work  for  this  research.  The  report  was  prepared  by 
Mr.  White  and  Dr.  Ebeling  under  the  direct  supervision  of  Mr.  Jones  and 
Mr.  Timothy  D.  Abies,  Acting  Director,  ITL. 

At  the  time  of  publication  of  this  report.  Dr.  James  R.  Houston  was  Director 
of  ERDC,  and  COL  John  W.  Morris  m,  EN,  was  Commander  and  Executive 
Director. 


The  contents  of  this  report  are  not  to  be  used for  advertising,  publication 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


V 


1  Introduction 


The  SOILSTRUCT-ALPHA  for  Personal  Computers  analysis  package 
comprises  three  separate  stages:  the  Visual  Modeler,  the  SOILSTRUCT-ALPHA 
Finite  Element  Analysis  Program,  and  the  Visual  Postprocessor.  This  report. 
Report  1  of  a  series,  describes  how  the  SOILSTRUCT-ALPHA  Visual  Modeler 
(i.e.,  a  preprocessor)  is  used,  applying  it  to  the  suggested  woik;  flow  of  a  typical 
project  as  an  example.  This  report  is  intended  to  serve  as  a  user’s  guide  and  not 
as  a  reference  manual  defining  each  specific  command.  Report  2  will  describe  the 
Finite  Element  Analysis  Program,  and  Report  3,  the  Visual  Postprocessor. 

1 .1  Overview  of  SOILSTRUCT-ALPHA  Finite 
Element  Program 

Engineering  application  of  soil-structure  interaction  modeling  requires  a 
balance  between  modeling  realism  and  simplicity.  Construction  procedures  have 
been  modeled  by  the  finite  element  method  for  over  25  years,  and  the  key 
ingredients  in  engineering  applications  are  well  known.  One  of  the  earliest 
successful  applications  of  soil-structure  interaction  analysis  was  performed  by 
Clough  and  Duncan  (1969;  Duncan  and  Clough  1971)  in  their  analysis  of  the  two 
reinforced  concrete  U-fi'ame  locks  at  Port  Allen  and  Old  River.  These  two  locks 
had  been  extensively  instrumented.  Prior  to  Clough  and  Dimcan’s  analyses,  the 
instrumentation  data  had  been  thought  to  be  unreliable  and  contrary  to  the 
perceived  understanding  of  the  behavior  of  locks  encountered  during  lock 
operation.  Clough  and  Duncan’s  study  showed  that  the  best  agreement  between 
results  computed  using  the  finite  element  method  and  those  obtained  through 
instrumentation  measurements  is  achieved  when  the  actual  construction  process  is 
simulated  as  closely  as  possible  in  the  analysis.  During  their  study,  Clough  and 
Duncan  developed  what  is  referred  to  as  the  backfill  placement  analysis  in  which 
the  loads  exerted  by  the  backfill  on  the  lock  wall  are  generated  automatically 
during  simulated  placement  of  backfill  behind  the  wall  (i.e.,  predetermined  earth 
pressure  force  distributions  between  the  soil  and  the  lock  are  not  specified).  This 
requires  that  the  soil  backfill  and  foundation  soil  strata  be  included  in  the  finite 
element  mesh.  This  procedure  involved  the  use  of  incremental  finite  element 
analysis  with  nonlinear,  stress-dependent,  stress-strain  behavior  for  the  soil.  An 
additional  requirement  is  that  interface  elements  be  incorporated  within  the  fimte 
element  mesh  to  allow  for  relative  movement  between  the  soil  and  structure. 
Clough  and  Duncan  developed  the  first  version  of  the  finite  element  program 
SOILSTRUCT  that  implements  this  engineering  modeling/analysis  philosophy. 
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SOILSTRUCT  has  been  used  successfully  in  the  past  decade  on  several 
engineering  projects  supported  by  field  observations.  Since  1969  several  versions 
of  SOILSTRUCT  have  been  developed  to  analyze  various  ^es  of  earth  retaining 
structures  or  to  analyze  specific  problems  that  were  not  envisioned  at  the  time  of 
Clough  and  Duncan’s  original  development.  One  of  these  updated  versions  is 
referred  to  as  SOILSTRUCT-ALPHA  and  is  the  subject  of  this  report. 

SOILSTRUCT-ALPHA  (Ebeling,  Duncan,  and  Clough  1990)  is  a  special- 
piupose  finite  element  program  for  two-dimensional  (2-D)  plane  strain  analysis  of 
soil-stmcture  interaction  problems.  SOILSTRUCT  calculates  displacements  and 
stresses  resulting  fi-om  incremental  construction,  backfilling,  excavation, 
dewatering,  rising  water  table,  and/or  load  application.  Nonlinear,  stress-path- 
dependent,  stress-strain  behavior  of  the  bacilli  was  approximated  in  the  finite 
element  analysis  using  the  tangent  modulus  method,  hi  the  tangent  modulus 
method,  new  values  of  tangent  moduli  are  assigned  to  each  soil  element  at  each 
increment  of  loading  (i.e.,  dewatering,  lock  construction,  and  backfilling)  or 
unloading  (i.e.,  excavation,  rising  water  table).  The  modulus  values  assigned  to 
each  element  are  adjusted  in  accordance  with  their  stresses  to  simulate  nonlinear 
behavior. 

Three  types  of  finite  elements  are  used  to  represent  the  behavior  of  different 
materials: 

a.  Two-dimensional  continua  elements.  A  2-D,  subparametric,  quadrilateral 
element  (QMS)  is  used  to  represent  the  soil  and  most  structural  materials. 
Structural  supports,  such  as  the  struts  or  tieback  components  of  an 
excavation  support  system,  are  typically  modeled  as  a  spring  support 
using  bar  elements.  However,  2-D  elements  have  been  used  to  model 
these  supports. 

b.  Interface  elements.  SOILSTRUCT-ALPHA  has  the  ability  to  model  file 
interface  region  between  the  soil  backfill  and  the  structure  using  interface 
elements.  This  important  feature  allows  for  the  movement  of  the  softer 
continua  elements  used  in  modeling  the  backfill  relative  to  the  movement 
of  the  stiffer  continua  elements  used  in  modeling  the  structure.  This 
element  is  defined  by  four  nodes,  with  each  of  the  two  pairs  of  nodes 
having  the  same  coordinates;  thus  this  element  has  no  thickness. 

c.  One-dimensional  bar  elements.  One-dimensional,  two-node,  bar  or 
spring  elements  are  used  to  model  the  behavior  of  a  variety  of  structural 
systems.  This  includes  the  modeling  of  structural  supports  such  as  braces 
or  tiebacks  or  the  modeling  of  reinforcement  placed  within  a  soil  backfill 

SOILSTRUCT  was  expanded  during  the  U.S.  Army  Corps  of  Engineers’  first 
Repair,  Evaluation,  Maintenance,  and  Rehabilitation  Research  Program  to  model 
the  loss  of  contact  between  the  base  of  a  wall  (a  lock  in  this  case)  and  its  rock 
foxmdation  using  a  procedure  called  the  ALPHA  method  (Ebeling,  Duncan,  and 
Clough  1990;  Ebeling  et  al.  1992).  The  ALPHA  method  was  extended  to  soil 
elements  by  Regalado,  Duncan,  and  Clough  (1992)  to  reduce  numerical 
inaccuracies  in  soil  elements  that  are  at  or  near  failme. 
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The  continua  elements  used  to  model  the  soil  and  the  soil-to-structure 
interface  elements  that  may  have  failed  in  shear  at  one  stage  of  loading  have  the 
ability  to  recover  their  shear  stiffiiess  and  shearing  resistance  as  a  result  of  an 
increase  in  confining  pressures  at  some  later  stage  of  loading  in  tins  version  of 
SOILSTRUCT-ALPHA.  Several  other  improvements  have  been  made  to  the 
material  models,  including  the  new  Gomez,  Filz,  and  Ebeling  (2000a,  2000b) 
extended  load/unload/reload  hyperbolic  model  for  interfaces,  and  to  the  numerical 
procedures  implemented  within  SOE^STRUCT-ALPHA  based  on  experience 
gained  at  the  U.S.  Army  Engineer  Research  and  Development  Center  in  con¬ 
ducting  soil-structure  interaction  analyses  of  different  types  of  Corps  structures. 

The  SOn.STRUCT-ALPHA  for  Personal  Computers  analysis  package 
comprises  three  separate  stages:  the  Visual  Modeler,  the  SOILSTRUCT-ALPHA 
Finite  Element  Analysis  Program,  and  the  Visual  Postprocessor.  Figure  1-1 
illustrates  the  relationships  among  the  separate  stages  of  the  SOILSTRUCT- 
ALPHA  analysis  package. 


1.1.1  Visual  Modeler 

The  Visual  Modeler  is  the  preprocessing  stage  for  SOILSTRUCT-ALPHA  for 
Personal  Computers.  The  Visual  Modeler  itself  comprises  two  programs,  the 
Mesh  Generator  (MeshGeri)  and  the  Load  Generator  (LoadGen).  The  Mesh 
Generator  is  a  computer-aided  design  program  for  user-controlled,  expedient 
generation  of  2-D  finite  element  meshes.  The  Mesh  Generator  is  a  line-based 
system  that  generates  three-  and  four-node  elements  automatically  with  nodes 
computed  fiom  the  intersection  of  line  segments.  To  expedite  the  generation  of 
mesh  lines,  the  user  defines  the  boundaries  and  material  areas  of  the  model  with 
precise  coordinates,  and  then  meshes  the  resulting  boundaries  with  tools  that 
allow  connecting  the  boundary  lines  with  an  arbitrary  number  of  lines  spaced  with 
a  user  control  based  on  Bezier-spline  methods.  The  Mesh  Generator  stores  its 
data  in  a  line-based  file  with  the  extension  .lin,  and  generates  an  element  file, 
designated  by  the  extension  .msh,  which  is  used  by  the  Load  Generator.  The 
Mesh  Generator  was  separated  fi'om  the  Load  Generator  because  the  Mesh 
Generator  was  a  more  general  program  and  could  be  used  for  other  2-D  finite 
element  modeling  work,  and  because  enhancements  could  be  added  more  easily, 
without  fear  of  impacting  the  load  addition  methods. 

The  Visual  Modeler’s  Load  Generator  is,  by  nature,  specific  to  the 
SOILSTRUCT-ALPHA  for  Personal  Computers  analysis  package.  The  Load 
Generator  is  designed  to  allow  the  user  to  quickly  and  easily  select  nodes  or 
elements  that  will  be  loaded  with  any  of  SOILSTRUCT-ALPHA’s  ten  load  types, 
and  to  accurately  place  bar  elements.  Plotting  the  elements  at  an  80  percent  scale 
enables  the  selection  of  nodal  points  relative  to  elements.  Selection  is 
accomplished  using  Region  Assignment,  where  the  user  draws  a  polygon  that 
surroimds  the  element  centers  or  nodal  points,  and  then  selects  file  attributes 
assigned  to  that  region.  Load  cases  are  viewed  with  their  appropriate  loads 
separately,  and  file  load  cases  can  be  added  or  removed  quickly  and  easily.  The 
Load  Generator  reads  in  the  .msh  file  exported  fi’om  the  Mesh  Generator  and 
stores  the  loads  in  an  .Ids  file.  The  program  also  generates  a  .dat  file  that  is  the 
input  for  the  SOILSTRUCT-ALPHA  Finite  Element  Analysis  Program. 
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Figure  1-1 .  Programs  for  SOILSTRUCT-ALPHA  and  their  associated  fiie 
extensions 
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1.1.2  SOILSTRUCT-ALPHA  Finite  Element  Program 

The  SOILSTRUCT-ALPHA  Finite  Element  Analysis  Program  is  a  stand¬ 
alone  executable  version.  It  requires  only  a  data  input  file  with  the  extension  .dat, 
which  is  output  fi'om  the  Visual  Modeler’s  load  generation  program.  It  also 
outputs  a  .pst  file,  which  will  be  used  by  the  Visual  Postprocessor. 


1.1.3  Visual  Postprocessor 

The  Visual  Postprocessor  is  used  to  view  the  results  stored  in  a  .pst  file  (i.e., 
element  stresses  and  nodal  point  displacements)  fi-om  the  nonlinear  finite  element 
analysis  of  the  soil-structure  interaction  problem  modeled  using  SOILSTRUCT- 
ALPHA  for  Personal  Computers.  This  visual  postprocessor  allows  the  user  to 
interpret  computed  results  visually. 


1.2  Modeling  Guidelines 

A  Visual  Modeler  user  may  find  the  following  guidance  helpfiil: 

a.  SOILSTRUCT-ALPHA  assumes  a  right-hand  coordinate  system.  It  is 
advisable  that  the  users  prepare  a  sketch  of  their  structure  in  this 
coordinate  system  prior  to  execution  of  the  Visual  Modeler. 

b.  In  a  SOILSTRUCT-ALPHA  construction  sequence  analysis  using  the 
backfill  placement  metiiod  of  analysis,  it  is  advisable  that  the  user  follow 
the  actud  geomehy  and  construction  sequence  as  closely  as  possible 
when  assembling  the  finite  element  model  and  in  the  subsequent 
incremental  analyses. 

c.  In  the  backfill  placement  analysis,  the  backfill  is  included  in  the  finite 
element  mesh,  and  the  forces  acting  on  the  wall  are  developed  through  a 
backfill  operation  simulation. 

d.  In  a  backfill  placement  analysis,  the  backfill  operation  is  simulated  by 
placement  of  the  backfill  in  a  series  of  lifts.  For  each  lift,  the  soil  placed 
generates  a  set  of  forces  that  are  applied  to  the  wall  and  previously  placed 
backfill  system.  With  an  equivalent  linear  stress-strain  soil  model  used, 
the  moduli  within  the  existing  soil  backfiU  are  adjusted  to  be  consistent 
with  die  current  level  of  loading.  In  a  backfill  placement  analysis,  two 
iterations  are  typically  specified  for  each  lift;  Ae  ALPHA  method  is  not 
implemented  until  the  second  computation  is  completed  per  lift. 

e.  Interfiice  elements  are  placed  in  the  following  situations: 

(1)  Between  the  backfill  and  the  structure  to  allow  for  relative  movement 
between  the  soil  and  the  structure. 

(2)  Between  the  backfill  and  the  rock  foimdation  (when  present)  to  allow 
for  relative  movement  between  the  soil  and  rock. 
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(3)  Between  the  structure  and  its  foundation,  be  it  a  soil  or  rock 

foundation,  to  allow  for  relative  movement  between  the  structure  and 
its  foundation. 

f.  The  ALPHA  method  was  developed  by  Ebeling,  Duncan,  and  Clough 
(1990)  using  data  computed  within  structure-to-rock  interface  elements  in 
a  numerical  algorithm  that  allows  for  base  separation  to  occur  between  the 
base  of  a  massive  concrete  lock  wall  and  its  rock  foundation.  Special 
requirements  are  placed  on  these  particular  interface  elements: 

(1)  The  base  of  the  lock  wall-to-rock  foimdation  interface  must  be 
horizontal. 

(2)  The  interface  elements  are  numbered  consecutively  starting  with  the 
interface  element  at  the  toe  of  the  wall  and  ending  with  the  interface 
element  at  the  heel  of  the  wall. 

Further  details  on  the  ALPHA  method  are  given  on  pages  64  through 
70  in  Ebeling  et  al.  (1992). 

g.  Cook,  Malkus,  and  Plesha  (1989,  page  578)  note  that  a  2-D  continua 
element  performs  best  if  its  shape  is  compact  and  regular.  They  also 
observe  that  an  element  tends  to  stiffen  and  lose  accuracy  as  its  aspect 
ratio  increases.  A  typical  recommendation  is  that  the  aspect  ratio  for  2-D 
continua  elements  be  restricted  to  5  to  1  for  meshes  analyzed  using  a 
linear  elastic  finite  element  method  of  analysis  to  compute  accurate 
results.  This  is  reasonable  guidance  for  linear  elastic  elements  that  are 
used  in  SOILSTRUCT-ALPHA  to  model  structures  made  of  concrete  or 
steel  and  rock  foundations.  However,  based  on  experience  by  the  second 
author  of  this  report,  this  general  guidance  is  too  generous  for  soil  regions 
of  interest  in  SOILSTRUCT-ALPHA  meshes.  For  example,  in  an 
incremental,  equivalent  linear  method  of  analysis,  it  is  recommended  that 
the  aspect  ratio  of  soil  elements  be  nearly  1  to  1  within  critical  regions 
such  as  at  the  soil-to-structure  interface  region  along  the  back  of  the 
retaining  structure.  Within  other  soil  regions  of  less  importance,  the 
aspect  ratio  can  be  increased.  Every  attempt  is  made  by  the  second  author 
of  this  report  to  limit  the  aspect  ratio  to  3  to  1  for  soil  elements. 

h.  Skewed  and  triangular  2-D  continua  elements  are  acceptable  for  linear 
elastic  elements  used  to  model  a  concrete  or  steel  structure  or  rock 
foundations.  (A  skewed  element  is  defined  as  the  comer  angle  between 
adjacent  element  faces  that  is  not  equal  to  90  degrees.  Note  that  at  every 
comer  of  both  square  rectangular  elements,  the  comer  angle  is  exactly  90 
degrees.)  Cook,  Malkus,  and  Plesha  (1989,  page  578)  observe  that  an 
element  tends  to  stiffen  and  lose  accuracy  as  the  comer  angles  become 
markedly  different  fl-om  one  another.  In  general,  the  more  skewed  a  2-D 
continua  element  is,  the  more  likely  its  stiffness  will  be  numerically  (that 
is,  artificially)  enhanced.  This  is  an  undesirable  effect.  Again,  an  element 
performs  best  if  its  shape  is  compact  and  regular.  This  general  advice 
holds  trae  for  meshes  analyzed  using  SOILSTRUCT-ALPHA. 
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i.  For  an  incremental,  equivalent  linear  method  of  analysis,  it  is 
recommended  that  the  use  of  triangular  elements  be  minimized  in  the 
mesh,  especially  within  soil  regions  of  interest.  Triangles  should  be 
avoided,  if  possible,  at  the  soil-to-structure  interface  region  along  the  back 
of  die  retaining  structure  and  other  critical  regions. 

j.  Clough  and  Duncan  (1969)  showed  that  the  accuracy  in  computed  results 
for  a  backfill  placement  analysis  improves  as  the  number  of  soil  lifts  used 
is  increased.  In  a  backfill  placement  analysis,  eight  soil  lifts  are  generally 
accepted  as  the  minimum  number  of  lifts  that  should  be  used  to  maintain 
sufficient  accuracy  in  computed  results.  Clough  and  Duncan  also  showed 
(in  their  Figure  25  for  “average  stresses”)  that  the  magnitude  of  the  error 
dso  depends  on  the  magnitude  of  the  hyperbolic  soil  parameter  exponent 
n,  with  larger  errors  associated  with  larger  values  of  n.  In  their  parametric 
evaluation,  the  values  for  n  were  between  0  and  1.1.  These  error  data 
were  generated  using  the  results  fi'om  SOILSTRUCT  analyses  and 
closed-form  solutions  for  deflections  of  a  one-dimensional  (1-D) 
representation  of  incremental  fill  placement. 

k.  Clough  and  Duncan  (1969)  also  showed  that  the  accuracy  in  computed 
results  for  an  excavation  analysis  improves  as  the  number  of  soil  layers 
used  is  increased.  They  developed  a  relationship,  given  in  their  Figure  24 
(for  “average  stresses”),  between  the  error  in  simulation  of  small  layer 
excavation  as  a  function  of  the  number  of  layers  used  in  the  analysis.  The 
magnitude  of  the  error  in  excavation  simulation  also  depends  on  the 
magnitude  of  the  hyperbolic  soil  parameter  exponent  n,  with  larger  errors 
associated  with  larger  values  for  n.  Their  Figure  24  was  generated  using 
the  results  from  SOILSTRUCT  analyses  and  closed-form  solutions  for 
deflections  of  a  1-D  representation  of  incremental  excavation.  No  general 
guidance  has  evolved  over  the  years  as  to  the  munber  of  incremental 
excavation  layers  that  should  be  used  to  maintain  sufficient  accuracy  in 
computed  results.  However,  the  second  author  of  this  report  attempts  to 
use  a  minimum  of  10  to  15  layers  in  an  excavation  analysis,  depending 
upon  the  structure  being  analyzed. 

1.3  Hyperbolic  Soil  Modeling  Parameter 
Guidelines 

SOILSTRUCT-ALPHA  uses  a  hyperbolic  representation  for  nonlinear  stress- 
strain  relationships  for  soil,  as  described  by  Duncan  and  Chang  (1970).  The  use  of 
this  mathematical  relationship  to  represent  the  primary  loading  curve  of  soil 
allows  for  the  development  of  a  mathematical  expression  for  &e  tangent  modulus. 
Recall  that  the  tangent  modulus  method  combined  with  the  incremental 
construction  method  of  analysis  accounts  for  the  nonlinear  stress-strain  behavior 
of  soil.  The  value  of  the  tangent  modulus  is  a  function  of  (a)  the  stresses  computed 
within  each  element,  (b)  the  Mohr-Coulomb  soil  strength  parameters,  and  (c)  a 
group  of  hyperbolic  parameters  that  are  obtained  from  triaxial  soil  test  data.  The 
stresses  computed  within  each  soil  element  are  characterized  in  terms  of  the 
confining  pressure  and  the  level  of  shearing  applied  to  the  element  of  soil.  The 
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confining  pressure  is  expressed  in  terms  of  the  minor  principal  stress  in  the 
hyperbolic  model.  The  level  of  shearing  applied  to  the  element  of  soil  is  expressed 
in  terms  of  the  principal  stress  difference.  The  values  for  confining  pressure  and 
the  level  of  shearing  applied  to  a  soil  element  will  change  during  each  iteration 
and  in  every  load  step  during  the  course  of  the  analysis.  Consequently,  the 
modulus  of  each  soil  element  is  adjusted  during  each  iteration  of  every  load  step 
during  the  course  of  the  analyses  in  accordance  with  the  stresses  in  the  element. 

A  Visual  Modeler  user  may  find  the  following  guidance  helpful: 

a.  The  hyperbolic  representation  for  nonlinear  stress-strain  relationships  for 
soil  requires  that  shear  test  data  be  obtained  fi'om  triaxial  tests  on  high- 
quality  imdisturbed  soil  specimens  recovered  fi-om  each  soil  stratum  in  the 
field.  Appendix  A  in  CJomez,  Filz,  and  Ebeling  (2000b)  outlines  the  steps 
involved  in  the  reduction  of  these  data  to  obtain  the  Mohr-Coulomb  soil 
strength  parameters  and  the  group  of  hyperbolic  parameters  fi’om  these 
triaxial  soil  test  data.  The  determination  of  the  Mohr-Coulomb  and 
hyperbolic  model  parameters  is  shown  for  isotropically  consolidated 
triaxial  tests  on  dense  Blacksburg  Sand  sheared  in  a  drained  mode  (i.e., 
isotropically  consolidated  drained  (ICD)  tests). 

b.  Over  the  years,  a  large  number  of  triaxial  shear  tests  for  different  soils 
have  been  performed  in  which  hyperbolic  stress-strain  model  parameters 
have  been  extracted.  A  number  of  these  tabulated  data  summaries  have 
been  collected  in  Appendix  A  of  Gomez,  Filz,  and  Ebeling  (2000b).  The 
tabulated  hyperbolic  parameter  values  data  contained  in  this  appendix  are 
not  intended  to  replace  site-specific  testing  on  high-quality  soil 
specimens,  but  are  for  use  as  interim  guidance  prior  to  availability  of 
laboratory  test  results.  These  data  can  be  used  as  an  aid  for  selecting 
hyperbolic  parameter  values  for  analysis  when  test  data  are  not  yet 
available.  It  has  been  observed  that  the  hyperbolic  parameters  correlate 
according  to  soil  type  and  with  standard  engineering  indices. 


1.4  Interface  Modeling  Parameters  Guidelines 

SOILSTRUCT-ALPHA  contains  the  Gomez,  Filz,  and  Ebeling  (2000a, 
2000b)  extended  load/unload/reload  hyperbolic  model  for  interfaces.  The  model  is 
based  on  the  Clough  and  Duncan  (1969)  hyperbolic  representation  for  nonlinear 
shear  stress-shear  deformation  relationship  for  primary  loading,  which  was 
extended  to  model  a  variety  of  stress  paths.  The  use  of  this  mathematical 
relationship  to  represent  the  primary  loading  curve  of  soil-to-structure  interfaces 
allows  for  the  development  of  a  mathematical  expression  for  the  tangent  interface 
stiffiiess.  Recall  that  the  tangent  stiffiiess  method  combined  with  the  incremental 
construction  method  of  analysis  accounts  fi)r  the  nonlinear  shear  stress-shear 
deformation  behavior  of  soil-to-structure  interfaces.  The  value  of  the  tangent 
interface  stiffiiess  is  a  fimction  of  (a)  the  normal  and  shear  stresses  computed 
within  each  interface  element,  (b)  the  Mohr-Coulomb  soil-to-structure  interface 
strength  parameters,  and  (c)  a  group  of  hyperbolic  parameters  that  are  obtained 
fi’om  soil-to-structure  interface  test  data.  The  stresses  computed  within  each 


8 


Chapter  1  Introduction 


interface  element  are  characterized  in  terms  of  the  normal  pressure  and  the  level 
of  shearing  applied  to  the  interface  element.  The  values  for  normal  pressure  and 
the  level  of  shearing  applied  to  an  interface  element  will  change  during  each 
iteration  and  in  every  load  step  during  the  course  of  the  analysis.  Consequently, 
the  stifiEuess  of  each  interface  element  is  adjusted  during  each  iteration  of  every 
load  step  during  the  course  of  the  analyses  in  accordance  with  the  stresses  in  the 
interface  element.  Note  that  the  Gomez,  Filz,  and  Ebeling  extended  load/unload/ 
reload  interface  model  does  not  require  any  material  property  values  in  addition  to 
those  considered  in  the  Clough  and  Duncan  (1969)  hyperbolic  interface  model. 

A  Visual  Modeler  user  may  find  the  following  guidance  helpful: 

a.  The  hyperbolic  representation  for  nonlinear  stress-strain  relationships  for 
soil  requires  that  shear  test  data  be  obtained  fi-om  stracture-to-soil 
interface  tests.  Appendix  B  in  Gomez,  Filz,  and  Ebeling  (2000b)  outlines 
the  steps  involved  in  tiie  reduction  of  these  data  to  obtain  the  Mohr- 
Coulomb  concrete-to-soil  interface  strength  parameters  and  the  group  of 
hyperbolic  parameters  fi’om  these  soil-to-structure  interface  test  data.  The 
determination  of  the  Mohr-Coulomb  and  hyperbolic  interface  model 
parameters  are  shown  for  concrete-to-soil  interface  tests  on  dense 
Blacksburg  Sand.  The  data  given  in  this  appendix  are  finm  interface  tests 
performed  in  the  Large  Displacement  Shear  Box  (LDSB)  at  Virginia 
Polytechnic  Institute  and  State  University,  Blacksburg.  The  LDSB  soil-to- 
concrete  interfaces  can  be  as  large  as  71 1  by  406  mm  (28  by  16  in.) 
allowing  a  maximum  interface  displacement  as  large  as  305  mm  (12  in.) 
(Shallenberger  and  Filz  1996).  Gomez,  Filz,  and  Ebeling  (2000a)  observe 
Aat  this  size  direct  shear  box  (DSB)  overcomes  two  of  the  limitations  of 
traditional  DSB:  The  maximum  relative  displacement  that  can  be  attained 
in  a  conventional  DSB  is  limited;  hence,  the  determination  of  the 
interface  residual  strength  becomes  difficult.  In  addition,  end  effects, 
induced  by  the  presence  of  the  rigid  walls  of  the  soil  container,  may 
induce  errors  in  the  results. 

b.  Over  file  years,  a  limited  number  of  experimental  studies  have  been 
conducted  to  obtain  data  on  the  behavior  of  soil-to-structure  interfaces  for 
different  soils:  Potyondy  (1961),  Clough  and  Duncan  (1969),  Peterson  et 
al.  (1976).  Their  tests  were  performed  using  a  slightly  modified,  standard 
size  DSB  in  which  a  concrete  specimen  occupied  one  of  the  halves  of  the 
shear  box.  Potyondy  (1961)  provides  Mohr-Coulomb  interface  shear 
strength  parameters  for  a  variety  of  structural  materials  and  soils.  Clough 
and  Duncan  (1969)  and  Peterson  et  al.  (1976)  provide  Mohr-Coulomb 
interface  shear  strength  parameters  for  the  interface  between  concrete  and 
different  types  of  soils  as  well  as  the  corresponding  hyperbolic  interface 
model  parameters  for  primary  loading.  Gomez,  Filz,  and  Ebeling  (2000a, 
2000b)  report  on  LDSB  interface  tests  between  concrete  and  a  variety  of 
granular  soils.  Table  3-8  in  Gomez,  Filz,  and  Ebeling  (2000b)  contains  a 
summary  of  hyperbolic  parameter  values  for  these  interface  tests.  The 
citations  listed  in  this  paragraph  for  tabulated  hyperbolic  parameter  values 
are  not  intended  to  replace  site-specific  testing,  but  are  for  use  as  interim 
guidance  prior  to  availability  of  laboratoiy  test  results.  These  data  can  be 
used  as  an  aid  for  selecting  hyperbolic  parameter  values  for  analysis  when 
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test  data  are  not  yet  available.  It  has  been  observed  that  the  hyperbolic 
interface  parameters  correlate  according  to  soil  type  with  standard 
engineering  indices  (e.g.,  Peterson,  et  al.  1976),  and  with  hyperbolic 
stress-strain  soil  parameters  (Table  3-9  in  Gomez,  Filz,  and  Ebeling 
2000b). 


1.5  Example  Structure 

An  example  of  a  complete  soil-structure  interaction  analysis  of  a  lock  wall 
retaining  earth  using  SOILSTRUCT-ALPHA  is  given  in  Ebeling  and  Wahl 
(1997).  Figure  1-2  shows  a  portion  of  the  cross  section  of  the  new  McAlpine 
roller-compacted  concrete  lock  wall  that  was  analyzed  using  SOILSTRUCT- 
ALPHA.  This  lock  structure  and  its  backfill  will  serve  as  an  example  for  the  step- 
by-step  assembly  of  a  complete  SOILSTRUCT-ALPHA  input  data  file  using  the 
Visual  Modeler.  The  finite  element  mesh  data  will  be  assembled  using  MeshGen, 
and  the  load  cases  to  be  applied  to  the  structure  are  assembled  using  LoadGen. 
Note  that  due  to  figure  size  limitations  for  producing  this  report,  the  rock 
foundation  that  was  included  in  the  origin^  analysis  is  excluded  in  this 
demonstration. 


Figure  1-2.  McAlpine  Lock  and  Dam  analysis  (Ebeiing  and  Wahl  1997) 
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1.6  Nomenclature 

For  the  purposes  of  this  report,  the  following  nomenclature  will  be  used 
Program  names  will  always  be  italicized. 

Definitions:  Definitions  will  be  surrounded  by  a  double-lined 
boimding  box.  _ 

Filename  extensions  will  be  preceded  by  a  dot  and  use  a  bold  typeface 
(i.e.,  .pst). 

Names  of  tools  will  be  in  bold. 
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2  MeshGen,  the  Mesh 
Generator 


As  stated  in  Chapter  1,  the  Mesh  Generator  is  a  computer-aided  design 
program  for  user-controlled,  expedient  generation  of  2-D  finite  element  meshes. 
This  program  is  a  line-based  system  that  generates  three-  and  four-node  elements 
automatically  with  nodes  computed  fi-om  the  intersection  of  line  segments.  To 
expedite  the  generation  of  mesh  lines,  the  user  defines  the  boundaries  and  material 
areas  of  the  model  with  precise  coordinates,  and  then  meshes  the  resulting 
boundaries  with  tools  that  allow  connecting  the  boundary  lines  with  an  arbitrary 
number  of  lines  spaced  with  a  user  control  based  on  Bezier-spline  methods.  The 
Mesh  Generator  stores  its  data  in  a  line-based  file  with  the  extension  .lin ,  and 
generates  an  element  file,  designated  by  the  extension  .msh,  which  is  used  by  the 
Load  Generator.  The  Mesh  Generator  was  separated  fi’om  the  Load  Generator 
because  the  Mesh  Generator  was  a  more  general  program  and  could  be  used  for 
other  2-D  finite  element  work,  and  because  enhancements  could  be  added  more 
easily  without  fear  of  impacting  the  load  addition  methods. 


Zooming:  Zooming  can  be  performed  at  any  time  by  click-dragging  a 
zoom  region  with  the  right  mouse  button.  To  return  the  zoom  to  view  the 
full  extent  of  the  model,  double-click  with  the  right  mouse  button. 
Zooming  comes  in  handy  when  hying  to  be  precise  with  that  assignment 
region,  and  especially  for  finding  that  closing  handle. 


2.1  Step  1 :  Initialize  the  Editor  for  the  Model 

In  this  step  the  user  needs  to  perform  the  following  tasks  to  set  standard  units 
and  extents  for  the  newly  created  model: 

•  Set  the  units  for  the  current  problem,  either  SI  or  non-SI. 

•  Set  the  drawing  limits  for  the  current  problem. 

•  Specify  the  snap  tolerance  of  the  units. 
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2.2  Step  2:  Define  the  Boundaries  of  the  Model 

Using  the  Line  and  Plumb  Line  tools  of  the  MeshGen  editor,  add  the 
boundary  lines  for  the  model.  These  boundary  lines  will  separate  different 
material  areas  and  the  complete  boundary  of  the  model. 


Line 


Plumb  Line 


2.3  Step  3:  Break  the  Material  Boundaries  into 
Four-  or  Three-sided  Regions 


After  the  boundary  definition  stage,  the  model  has  many  different  material 
areas,  some  of  which  might  be  complex.  It  is  best  to  take  these  regions  and  break 
them  into  manageable  parts.  Use  the  Line,  Break  Line,  Weld  Line,  Fuse  Line, 
Extend  Line,  and  Plumb  Line  tools  to  simplify  the  regions  to  help  with  the 
meshing  tools. 


Line 

Break  Line 
Weld  Line 
Fuse  Line 
Extend  Line 
Plumb  Line 


Chapter  2  MeshGen,  the  Mesh  Generator 


13 


Once  everything  is  in  place,  use  the  menu  item  Split  All  Lines  from  the  Edit 
menu  to  make  line  segments  that  can  be  conveniently  used  in  Step  4  (Mesh  the 
Regions).  This  command  separates  continuous  lines  into  line  segments  at  all  the 
intersection  points  where  lines  cross  one  another. 


Edit 
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2.4  Step  4:  Mesh  the  Regions 


Using  the  Combine  Lines  tool,  select  two  lines  to  mesh  between.  Assign  the 
number  of  lines  and  the  method  (linear  or  nonlinear)  used  to  connect  the  selected 
lines.  Continue  until  the  entire  region  is  meshed. 


Combine  Lines 
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2.5  Step  5:  Assign  Interface  Lines 


Certain  lines  can  be  marked  to  become  interface  elements  when  the  mesh  is 
generated.  These  lines  can  be  marked  by  selecting  them  (using  the  left  mouse 
button  and  dragging  a  region  that  overlaps  or  surrounds  the  line)  and  by  choosing 
Set  Selected  Lines  as  Interface  Elements  from  the  Edit  menu.  There  is  also  a 
method  to  reset  a  line. 


Edit 

l 


File  DMifigure  Edit  ^ely2e 
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2.6  Step  6:  Generate  Mesh  (which  will  be  input 
into  the  Load  Generator) 

From  the  Analyze  menu  of  MeshGen,  choose  Generate  Mesh.  The  mesh 
generation  automatically  builds  nodes  and  elements  from  the  line  drawings  of 
MeshGen. 


Anajyze 


File  Configure  Edl  ^ajyze 


2.7  step  7:  Optimize  Mesh  (optional) 

Perhaps  the  project  takes  a  long  time  to  run  under  SOILSTRUCT-ALPHA,  or 
even  does  not  run  because  of  memory  constraints.  Sometimes  the  best  solution  is 
to  renumber  the  nodes  in  the  mesh.  An  efficient  node  numbering  optimizer  can 
reduce  the  size  of  the  matrices  involved  in  finding  a  solution.  Because  it  is  faster 
to  operate  on  smaller  matrices,  the  newly  renumbered  mesh  is  more  time  and 
space  efficient.  However,  the  new  node  configuration  is  usually  nonintuitive, 
forcing  the  user  to  search  for  the  node  he  or  she  would  like  information  on.  For 
this  reason,  there  are  three  ways  to  number  nodes  for  bandwidth  minimization  in 
MeshGen. 

Under  the  Analyze  menu  of  the  MeshGen  program,  there  is  an  option  labeled 
Advancing  Wavefront  Optimization.  Of  the  renumbering  methods,  this  option 
provides  the  most  efficient  and  least  intuitive  renumbering  for  nodes.  This 
renumbering  scheme  takes  connectivity  into  account  when  it  tries  to  optimize  the 
mesh. 

The  other  two  options  deal  only  with  sorting  the  nodes,  whether  in  the  X-  or 
Y-direction  first,  and  in  the  complementary  direction  second.  If  a  mesh  is  taller 
than  it  is  wide,  it  is  usually  a  good  idea  to  sort  in  the  X-direction  first.  If  the 
model  is  wider  than  it  is  tall,  then  it  will  generally  be  better  to  sort  in  the 
Y-direction  first. 
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2.8  Step  8:  Verify  Mesh  -  Node  Numbers  and 
Element  Numbers 


From  the  Edit  menu,  choose  Plot  Mesh  to  see  the  newly  generated  mesh. 
Selecting  Element  Numbers  or  Node  Numbers  plots  the  element  and  node 
numbers,  respectively.  If  the  elements  and  nodes  are  not  in  the  order  expected,  go 
to  the  Configure  menu  item  Sort  Mode  to  look  at  different  options.  If  a  different 
option  is  selected,  the  mesh  will  have  to  be  generated  again.  To  aid  in  the 
process,  using  the  right  mouse  button  to  drag  an  area  will  cause  the  screen  to 
zoom  on  that  region.  Double-clicking  the  mouse  will  restore  the  full  view.  (This 
option  is  available  throughout  the  program.) 
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2.9  Step  9:  Refine  and  Reiterate 

Save  the  .lin  file  from  the  File  menu  and  modify  the  mesh  any  time.  With  the 
editing  tools,  lines  can  be  deleted,  created,  or  edited  over  and  over  again,  and  the 
mesh  regenerated. 

File 

i 


File  Configure  Edit  Analyze 


Select  Lines 
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LoadGen,  the  Load 
Generator 


As  stated  in  Chapter  1,  the  Visual  Modeler’s  Load  Generator  is,  by  nature, 
specific  to  the  SOILSTRUCT-ALPHA  for  Personal  Computers  analysis  package. 
The  Load  Generator  is  designed  to  allow  the  user  to  quickly  and  easily  select 
nodes  or  elements  that  will  be  loaded  with  any  of  SOILSTRUCT-ALPHA’s  ten 
load  types,  and  to  accurately  place  bar  elements.  Plotting  the  elements  at  an 
80  percent  scale  enables  the  selection  of  nodal  points  relative  to  elements. 
Selection  is  accomplished  using  Region  Assignment,  where  the  user  draws  a 
polygon  that  surrounds  the  element  centers  or  nodal  points,  and  then  selects  the 
attributes  assigned  to  that  region.  Load  cases  are  viewed  with  their  appropriate 
loads  separately,  and  the  load  cases  can  be  added  or  removed  quickly  and  easily. 
The  Load  Generator  reads  in  the  .msh  file  exported  from  the  Mesh  Generator  and 
stores  the  loads  in  an  .Ids  file.  The  program  also  generates  a  .dat  file  that  is  the 
input  for  the  SOILSTRUCT-ALPHA  Finite  Element  Analysis  Program. 


Assignment  regions:  Assignment  regions  are  created  by  left-clicking 
with  the  mouse  to  create  a  polygon  encompassing  either  element 
center(s)  or  nodal  point(s),  depending  on  whatever  the  load  condition 
affects.  The  assignment  region  must  be  closed  to  be  finished,  so  the 
first  point  selected  for  the  polygon  has  a  handle  (drawn  as  a  box) 
associated  with  it.  Left-clicking  with  the  mouse  inside  the  box  closes 
the  polygon  and  finishes  the  assignment  region.  Any  number  of  points 
can  be  used  to  create  an  assignment  region,  and  there  are  no  restrictions 
on  shape. 


Finishing  or  canceling  an  operation:  Some  operations  that  can  be 
performed  under  the  Load  Assignment  preprocessor  stage  do  not  have 
a  definite  ending  point,  and  some  operations  might  need  to  be 
cancelled  before  they  are  finished.  In  those  situations,  the  arrow  button 
on  the  toolbar  to  the  right  of  the  primary  display  will  display  the  word 
“Finished”  or  “Cancel”  and  the  operation  can  be  concluded  by  clicking 
on  the  button. 
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Zooming:  Zooming  can  be  performed  at  any  time  by  click-dragging  a 
zoom  region  with  the  right  mouse  button.  To  return  the  zoom  to  view 
the  full  extent  of  the  model,  double-click  with  the  right  mouse  button. 
Zooming  comes  in  handy  when  trying  to  be  precise  with  that 
assignment  region,  and  especially  for  finding  that  closing  handle. 


3.1  Step  1:  Load  the  Mesh  created  in  MeshGen 

From  the  File  menu,  choose  the  Open  (Input  Mesh)  option  and  load  the 
mesh  created  in  the  Mesh  Generation  program. 


File 
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3.2  Step  2:  Assign  Initial  Material  Properties  for 
2-D  and  Interface  Elements 


From  the  Material  Definition  menu,  choose  either  the  2D  Material  or 
Interface  Material  submenu.  The  dialog  that  appears  allows  for  the  entry  of  the 
specific  properties  of  each  material  that  will  be  used  in  the  model.  The  arrow 
buttons  next  to  the  material  number  allow  the  user  to  cycle  through  each  material 
and  set  its  properties  in  turn. 


Material  Definition 


2D  Material  Number:  [T-tj  Color:  j 

Material  Name:  |Sfl2D) 

Material  behavior  (IDRAIN):  I  Patch  Element  ^ 


Poisson's  ratio  before  failure  (GUI): 

Pijisson'i  i'btiij  [GL'r: 

Tnta.;  n f  b L- 1 ; y ^ t it  J ] ; ] . '  •.'..•e ■  c; r. r  ( P - ^.-i; . 

FciifiiP  ra'rn 

Friction  angle  in  degrees  (PHI): 

Modulus  exponent  (n)  (><XP): 

lifL'V  r  i.f i:--)  i;-  P.  ''FF''. 

Tangent  modulus  or  Young's  modulus  (E): 


Interface  Material  Number:  fn  -r^  Color:  | 

Interface  Material  Name:  [Air  [Interface) 


Interface  friction  angle  in  degrees  (PHI): 

Min.  value  for  interface  shear  stiffness  (RKS): 

Interface  cohesion  (COJ): 

Failure  ratio:  (FRJ) 

Tensile  strength  of  interface  (TENSJ): 

Adjacent  2D  element  material  type  number  (lADJMT); 
Interface  modulus  number  (RKJ): 

Interface  modulus  exponent  (XKPJ): 
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3.3  Step  3:  Assign  Fixities  for  the  Modei 

Fixities  are  assigned  to  nodal  points  using  assignment  regions.  From  the 
Load  Case  menu  choose  the  Assign  Fixities  menu  item  and  draw  out  an 
assignment  region  encompassing  the  nodal  points  to  be  fixed.  After  the 
assignment  region  has  been  selected,  a  dialog  will  appear  prompting  for  the  type 
of  fixity  to  apply. 


Load  Case 

^ . ^ . . i  . . .  . . .  . 

Re  Configuie  Edit  Material  Definition  Load  Case  Computation  Parameteis 
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Current 
Load 
Case  : 


3.4  Step  4:  Assign  Initial  Materials  to  the  Mesh 

The  initial  materials  for  the  mesh’s  2-D  and  interface  elements  are  assigned 
using  assignment  regions.  From  the  Load  Case  menu,  select  either  Assign  2D 
Element  Material  or  Assign  Interface  Materials,  encompass  an  assignment 
region,  and  select  a  material  from  the  dialog. 

Jpri i i  'i  ?  i » I  * •'  f ;  M :  ;■  1  li  i  L f  iit-i  - i  i .  .  ■  ■ , .  ■  -Tnt  x| 

Re  Configure  idt  Meterial  Definition  Load  Case  Computation  Parameters  1 


Current  ■ 
Load  i 
Case 
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3.5  Step  5:  Create  Initial  Bar  Elements  (Optional) 


Bar  elements  can  be  added  to  the  mesh  by  selecting  the  Bar  Element  button 
from  the  toolbar  to  the  right  of  the  primary  viewing  screen.  Each  time  the  button 
is  pressed,  one  bar  element  can  be  added  to  the  mesh.  Since  bar  elements  are 
connected  to  a  certain  node  of  an  element,  the  node  must  be  selected  by  click¬ 
dragging  with  the  left  mouse  button  over  the  shrunken  element  plot’s  nodal 
position  (i.e.,  the  appropriate  comer  of  the  shmnken  element).  When  both  end 
points  are  selected,  a  dialog  pops  up  asking  for  attributes  of  the  bar  element. 


3.6  Step  6:  Create  Load  Cases 

Load  Cases  have  to  be  created  on  an  individual  basis.  Selecting  Add  Load 
Case  from  the  Load  Case  menu  creates  a  new  load  case  for  the  current  model. 
The  current  load  case  is  accessed/changed  through  the  Current  Load  Case  control 
on  the  toolbar  to  the  right  of  the  primary  viewing  screen.  Individual  load  cases 
can  also  be  deleted. 

Load  Case 


i 
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3.7  Step  7:  Add  Regions  and  Conditions  for  Each 
Load  Case  (Such  as  Fili  Regions,  Upiift 
Conditions,  etc.) 

By  using  the  Load  Case  menu,  up  to  three  types  of  load  conditions  can  be 
added  to  each  load  case  in  the  model.  This  is  a  SOILSTRUCT-ALPHA 
limitation,  but  it  is  strictly  enforced  in  LoadGen.  When  a  load  type  exists  in  a 
load  case,  the  load  type  menu  item  is  flagged  for  visual  confirmation.  Most  load 
types  can  be  entered  using  the  region  assignment  method,  with  a  few  exceptions 
(i.e.,  bar  elements  and  seepage). 


Load  Case 


•;  I,  ‘  y.-  ^  ’■  k  '  ‘ » 'V'-  -- 1  k  -  v <  bi  libi?  i’.f  M  k! ''  ^’^5  b  ft t  •  J  f  t’< <  ik  •  j  i.  t  if* ( J-i  T-iH  f p- 

Rle  Configure  Etfi  Maletial  Definition  Load  Case  Computation  Parameters 


3[M}x] 


Current 
Load  ; 
Case 

□Iffl 


Bar  elements  can  be  added  in  a  method  similar  to  creating  the  initial  bar 
elements.  Bar  elements  can  also  be  deleted  by  click-dragging  with  the  left  mouse 
button  to  select  the  bar  element  to  be  deleted. 

Seepage  is  defined  by  water  levels.  When  Define  Seepage  has  been  chosen 
from  the  Load  Case  menu,  a  blue  line  representing  water  levels  is  drawn  over  the 
model,  which  has  been  faded  to  give  a  visual  cue  that  the  program  is  in  Define 
Seepage  mode.  Left  mouse  clicking  brings  up  a  dialog  box  with  the  current 
mouse  position  already  entered.  If  the  left  mouse  has  been  click-dragged  to 
encompass  a  nodal  point,  then  that  nodal  point  is  entered  into  the  dialog.  When 
the  dialog  is  accepted,  the  blue  line  is  changed  to  show  how  the  water  level  is 
changed.  If  the  points  entered  are  beyond  the  extents  of  the  model,  then  the  point 
is  ignored. 
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3.8  Step  8:  Create  Spreadsheet  Output  (Optional) 

SOILSTRUCT-ALPHA  has  an  option  to  allow  the  user  to  output  data  in  a 
format  that  spreadsheet  programs  can  read  easily.  This  option  allows  the  user  to 
quickly  manipulate  and  visualize  data  for  specific  nodes,  elements,  and  stresses 
along  edges  of  an  element. 

To  use  this  option  the  user  needs  to  select  which  data  he  or  she  desires  to 
output.  Selecting  Plot  Key  from  the  Spreadsheet  Output  submenu  of  the 
Computation  Parameters  menu  causes  a  dialog  to  appear.  This  dialog  has  a  list 
of  options  for  data  output,  either  stresses  for  elements,  displacements  for  nodes, 
normal  and  shear  stresses  for  interface  element  edges,  or  horizontal  and  vertical 
stresses  for  element  edges. 

After  a  plot  key  has  been  selected  for  output  data,  the  user  needs  to  select  the 
nodes  and  elements  for  output.  Selecting  Group  Edit  from  the  Spreadsheet 
Output  submenu  of  the  Computation  Parameters  menu  starts  the  node/element 
selection  mode.  The  user  can  then  create  an  assignment  region  about  the 
nodes/elements  that  he  or  she  wants  to  use. 

If  the  user  has  selected  stresses  for  2-D  elements,  then  any  2-D  element 
handles  in  the  assignment  region  output  those  elements.  Likewise,  if  the  user  has 
selected  displacements  for  nodes,  then  all  nodes  in  the  assignment  region  are 
output.  However,  if  the  user  selected  edge  information,  then  an  assignment  region 
must  surround  an  element  handle  and  two  nodes  along  one  edge  of  that  element 
before  that  information  is  output.  If  more  than  two  nodes  are  surrounded  by  that 
region,  then  each  set  of  nodes  for  each  edge  selected  is  output. 


Computation  Parameters 
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3.9  Step  9:  Set  Parameters  for  Analysis  (Optional) 

SOILSTRUCT-ALPHA  gives  the  user  the  power  to  affect  how  loading 
calculations  are  performed.  The  user  has  the  capability  to  change  parameters  for 
base  separation  and  parameters  used  to  check  excessive  incremental 
displacements  if  the  alpha  method  will  be  used. 

Alpha  parameters  for  base  separation  are  applied  by  selecting  the  Base 
Separation  Parameters  option  from  the  Computation  Parameters  menu. 

When  this  option  is  selected,  the  user  will  be  prompted  for  the  maximum  number 
of  iterations  using  the  alpha  method  for  each  load  step  and  for  the  type  of  stress 
increment  to  be  applied  using  the  alpha  method.  After  these  parameters  have 
been  entered,  the  user  will  be  required  to  identify  the  base  interface  element  group 
for  separation  using  the  alpha  method.  When  the  alpha  method  is  invoked,  the  ' 
interface  elements  forming  the  base  of  lock  wall-to-rock  foundation  interface  must 
be  horizontal.  The  LoadGen  program  will  automatically  form  interface  element 
groups  by  collecting  linear  interface  element  segments.  To  select  a  group,  a 
region  must  be  drawn  encompassing  the  group  (or  groups).  Selecting  the  base 
group  requires  the  user  to  select  the  single  horizontal  group  for  separation.  After 
the  user  has  selected  the  base  group,  he  or  she  will  be  prompted  to  select  other 
groups  for  which  resultant  forces  and  their  points  of  application  will  be  computed. 
Double-clicking  after  the  resultant  groups  are  selected  returns  the  user  to  the 
default  editing  mode. 

Control  parameters  for  checking  excessive  incremental  displacements  can  be 
applied  by  the  user  by  selecting  the  Excessive  Incremental  Displacement 
Parameters  option  from  the  Computation  Parameters  menu.  This  option  is 
used  when  the  alpha  method  is  invoked.  This  menu  selection  causes  a  dialog  box 
to  appear,  prompting  the  user  for  information  about  the  excessive  incremental 
parameters.  The  user  is  prompted  for  the  type  of  incremental  displacement  check 
to  perform,  the  initial  load  number  for  this  check,  tolerance  values  for  the  X-  and 
Y-directions,  and  whether  the  program  will  stop  when  a  tolerance  is  exceeded. 

After  the  dialog  box  has  been  completed,  the  user  will  be  allowed  to  select 
nodal  points  for  which  the  incremental  displacement  check  is  performed.  The 
user  creates  a  region  enclosing  the  nodes  to  check  or  double-clicks  to  choose  all  of 
the  nodes.  If  a  region  has  been  drawn,  all  nodal  points  in  that  region  are  checked. 


3.10  Step  10:  Set  Restart  Information  (Optional) 

SOILSTRUCT-ALPHA  provides  the  user  with  the  option  of  continuing  a 
previous  analysis.  In  order  for  the  user  to  use  this  option,  he  or  she  must  specify 
that  a  restart  is  taking  place  and  specify  the  analysis  results  restart  file  name  of  the 
analysis  to  be  restarted.  Access  the  Input  Restart  File  submenu  from  the  File 
menu  to  give  the  restart  information. 
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3.11  Step  11:  Save  the  Loads  into  an  .Ids  File 


When  the  load  cases  have  been  entered  for  a  specific  model,  the  load  cases 
can  be  saved  in  a  separate  file  for  easy  retrieval  and  editing.  From  the  File  menu, 
select  Save  (Loads).  The  load  cases  will  be  saved  with  the  same  name  as  the 
mesh  file,  but  with  the  extension  .Ids  .  Bear  in  mind  that  if  the  mesh  has  changed, 
the  load  case  file  associated  with  it  may  be  invalid.  However,  if  the  mesh  has  not 
changed,  then  its  load  cases  can  be  reopened  and  edited  easily  for  refinement  of 
the  analysis. 


File 

i 
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3.12  Step  12:  Create  SOILSTRUCT-ALPHA  Input 
File 

This  is  the  final  step  of  the  preprocessor  stage.  When  the  load  cases  have 
been  entered  for  a  specific  mesh,  and  the  output  options  selected,  a  valid 
SOILSTRUCT-ALPHA  input  file  can  be  created  for  processing.  Select  Create 
SOILSTRUCT-ALPHA  Input  from  the  File  menu  to  create  a  .dat  file  for  input 
into  the  SOILSTRUCT-ALPHA  processor. 
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Summary 


This  report  describes  how  the  SOILSTRUCT-ALPHA  visual  modeler  (i.e.,  a 
preprocessor)  is  used,  applying  it  to  a  suggested  work  flow  of  a  typical  project. 
The  SOILSTRUCT-ALPHA  for  Personal  Computers  analysis  package  comprises 
three  separate  stages:  the  Visual  Modeler,  the  SOILSTRUCT-ALPHA  Finite 
Element  Analysis  Program,  and  the  Visual  Postprocessor.  Figure  4-1  illustrates 
the  relationships  among  the  separate  stages  of  the  SOILSTRUCT-ALPHA 
analysis  package.  The  Visual  Modeler  is  intended  to  expedite  the  mesh  con¬ 
struction  and  load  generation  for  the  typical  SOILSTRUCT-ALPHA  input  file. 
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Figure  4-1.  SOILSTRUCT-ALPHA  Package  for  Personal  Computers 
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14.  (Concluded). 

The  Visual  Modeler  is  designed  to  facilitate  the  development  of  a  finite  element  mesh  for  SOILSTRUCT-ALPHA  in  a  timely 
manner.  The  modeler  comprises  two  parts,  a  computer-aided  design  tool  for  the  actual  creation  and  refinement  of  a  mesh,  including 
four-noded  two-dimensional  continua  and  interface  elements,  and  a  computer-aided  design  tool  for  the  application  of  loads  to  the 
mesh  and  the  creation  of  one-dimensional,  two-noded  bar  elements.  These  tools  expedite  the  creation  and  refinement  of  the  often 
extensive  input  data  to  SOILSTRUCT-ALPHA. 
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